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ABSTRACT 



Photoelastic methods using frozen stresses are very useful for 
experimentally determining stress distributions in bodies having 
complicated geometries for which no analytical solution is known. 

The method is applied to cylindrical pipe fittings having discontinui- 
ties and subjected to hydrostatic internal pressure loading. The equip- 
ment used, the method of model construction and the calibration of the 
model material are described. Two models having different geometries 
are tested and the results are compared with theory and with each other. 
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NOMENCLATURE 



English Letter Symbols 



a 


inside radius 


a 

m 

t 


mean radius 

* 

Angstrom units 


b 


outside radius 


E 


modulus of elasticity 


f 


photoelastic material fringe constant 
for stress 


h 


wall thickness or beam height 


I 


moment of inertia (second moment of area) 
about centroidal axis 


M 


moment 


N 


fringe order 


n 


total number of fringes present 


P 


force 


Pi 


internal pressure 


r 


radius, radial cylindrical co-ordinate 


t 


thickness of a slice or beam 


z 


.longitudinal distance, longitudinal cylindrical 
co-ordinate 



Greek Letter Symbols 



B 


parameter used in shall bending theory 

angular position, angular cylindrical co-ordinate 


h 


Poisson's ratio 


7T 


the constant 3. 14159. • • • • • • 


<r 


stress (general) 
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<r& 

<Tr 



general normal stresses 
circumferential normal stress 
radial normal stress 
longitudinal normal stress 
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Section 1. Introduction 



The use of photoelastic techniques to study three-dimensional models 
is not new; **stress freezing** has been known and applied for about thirty 
years. In recent years these methods have received renewed interest from 
many investigators because of the relative ease with which useful results 
can be obtained for problems with complicated geometries or having complex 
systems of loading, and for which analytical solutions are cumbersome or 
are not known. 

It was decided to investigate stress distributions in a cylindrical 
body of revolution containing discontinuities. The specific geometry 
chosen was that of a straight coupling joining two short sections of 
pipe. The loading selected was internal hydrostatic pressure. This 
problem was selected because it has rather wide engineering applications 
and because failures often occur near the joints in piping systems. 

Section 2. Model Design 

The exact geometry selected was that of a cylindrical, socket ended 
straight coupling joining two short sections of pipe. The pipe size 
corresponded to l-% inch ASTM Schedule 80 pipe having an inside diameter 
of 1.500 inch and an outside diameter of 1.900 inch. The matching 
coupling was 2-3/4 inches long, with an inside diameter of 1.500 inch 
and an outside diameter of 2.200 inches. The sockets were specified to 
be one-half inch deep, with an inside diameter appropriate to provide 
a slightly loose sliding fit with the pipe, or about 1.905 inches. The 
arrangement of the components is shown in figures 1 and 2. 

One model followed the above plan exactly; it will hereafter be 
referred to as Model II. The other model was modified by the addition 
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gure 1. Model I Assembly Drav/in; 
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Figure '2. Model II Assembly Dra.v/ing 



of a circumferential groove cut in the midsection of the coupling to a 
root diameter of 2.000 inches. This groove had a circular profile of 
1 inch radius. This modified model will be referred to in the text and 
figures as Model I. 

The pipe sections were jointed to the couplings by cementing to 
create a joint similar in appearance to silver brazed pipe joints. One 
end was modified slightly by the addition of a 45 degree bevel in order 
to simulate a comer filled with a non-flowing brazing alloy. 

The dimensions selected for the models were convenient to fit the 
oven and the polariscope field of view. An additional factor considered 
in sizing the pipe stubs was the presence of the end closures. These 
were almost rigid constraints on the pipe when the loading was applied, 
creating localized disturbances consisting of bending moments and 
shearing forces. According to the theory for bending of cylindrical 
shells as presented by Timoshenko and Woinowsky-Krieger, the disturb- 
ances die out with increasing distance along the cylinder. The theory 

X 7T 

states that effects are negligible for 2 . / — , where 

_ 3 a. A 
k" 

and where Poisson's ratio, a = mean radius of the cylinder 

/ ® 

h « wall thickness of the cylinder 

7T 

For the selected pipe dimensions, this value of z is ^ or 1.05 
inches. [5]^ Therefore, it was inferred that the end disturbances 
would not extend to the vicinity of the joint. The "safe” length cal- 
culated here can only be an estimate inasmuch as the theory was developed 
for axi- symmetric shells in which the original stress system was bi- 
axial, containing only membrane stresses, and in which the final stress 




Clumbers in square brackets refer to publications listed in the 
bibliography. 
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system created by the localized disturbances was composed of shell bend- 
ing stresses super-imposed upon the original system. The selected geo- 
metry has a triaxial stress system and other deviations from the assump- 
tions of the above theory; these are more fully described in Section 7. 

Section 3. Model Construction 

It was decided to machine the model components from commercially 
cast stock in order to avoid the difficulties associated with casting. 

The material selected for the model components was Hysol CP5-4290 cast 
epoxy in rod and tube form. The matching cement used for assembly of 
the components was Hysol R9-2039 liquid resin and H2-3475 hardener. 

These are stock items available from Hysol Corporation, 1100 Seneca 
Avenue, Olean, New York. The manufacturer's technical data for Hysol 
4290 is shown in Appendix I, page 46. 

The pipe sections were machined from two inch diameter tubing; little 
machining was required. The coupling, however, had to be made from three 
inch rod stock. This was a straightforward but lengthy process because 
epoxy must be machined using light cuts to avoid excessive heat genera- 
tion. (Figure 3 shows the rough stock, two pipe sections, two couplings 
and two calibration specimens;)) 

The pipe and coupling sections were joined using the liquid resin 
and hardener mixed in the proportion of two parts resin to one part 
hardener, by weight. It was at>plied liberally, squeezed to force air 
out, and the excess was wiped from the inner and outer surfaces of the 
assembly. Both end closures were formed by casting the liquid resin 
and hardener mixture directly into the open ends of the assembly, one 
end at a time, using aluminum foil molds. A short section of flared 
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copper tubing was molded into place in one end closure in order to 
provide a pressurizing connection for loading the model. A typical com- 
pleted assembly is shown in figure 4. 

All work with the cement required that the mating surfaces be 
thoroughly degreased. This cleaning was accomplished by dipping in dry 
cleaning solvent, drying in air, and finally wiping with acetone. After 
the surfaces had been prepared, the required quantity of cement was 
mixed, avoiding introduction of air bubbles as much as possible. Pot 
life after mixing was short--about fifteen minutes. The resin-hardener 
mixture reacted exothermically with a temperature rise dependent upon 
the quantity of mixture present in bulk. Thick sections heated them- 
selves enough to accelerate the curing process; curing was complete in 
about an hour. Thin sections, however, required an overnight cure at 
room temperature, or about one hour at 180®F. 

Section 4. Material Calibration 

Several methods of calibration may be employed in stress freezing. 
The method selected was to use a beam in pure bending as described in 
Handbook of Experimental Stress Analysis . [3] This type of specimen 
has a disadvantage of requiring fairly large pieces of stock, but it 
could be easily made, was easily loaded in available fixtures, could 
be tested in the same container with a model, and provided a number of 
data points with a single test. 

The length of the calibration beams was limited by the inside dia- 
meter of the oil container in which the beams were suspended for the 
stress freezing cycle. The beam height was selected to avoid the 
possibility of lateral buckling. Burke suggested a height of one-half 
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Inch. [1] This size was used for this project. 

The required specimens were fabricated by first bandsawing two 
slabs about one-third inch thick from an eight inch length of the three 
inch diameter epoxy rod. The slabs were machined with a fly cutter to 
obtain flat, parallel sides. Specimen blanks were cut slightly over- 
size and both surfaces were sanded by hand until thickness variation was 
within .001 inch. The holes were drilled, a template was attached, and 
the final machining of the edges was done with a high speed rotary file. 

The drawing of the calibration beam. Figure 5, shows the nominal 
dimensions. The exact dimensions of height and thickness were deter- 
mined for each individual specimen by measurement with a micrometer 
caliper and were tabulated for later use in the calibration calculations. 

Bach calibration specimen was subjected to a separate stress freez- 
ing cycle while it was immersed in mineral oil contained in a tank with- 
in the oven. The calibration specimen was either alone or it accompanied 
one of the models being tested. The reason for this arrangement was to 
Insure that the calibration specimen followed the same temperature pro- 
gram as the model in order to obtain a material fringe value best re- 
presenting the model material. 

The gravitational body force on a calibration beam is small. A Speci- 
mdn submerged in oil experiences an upward buoyant force which tends to 
reduce the effect of the body force still further, thus serving to re- 
duce the error Incurred by ignoring the additional bending moment caused 
by this force. A check calculation (Appendix I) shows that the error 
resulting from neglecting the apparent weight of the calibration beam 
is about one per cent. A similar effect was of course acting on the 
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Figure Calibration Specimen Drawing 



calibration loading weights; this was taken Into account simply by weigh- 
ing them while they were Immersed in hot oil to closely approximate the 
test conditions. 

A scale factor for fringe measurements on the photographs of the 
fringe patterns of the calibration specimens was established Individual- 
ly for each beam. Reference marks were scribed on the side to be photo- 
graphed and the separation distance was measured with a tool-maker's 
microscope. Measurement of the corresponding distance on the fringe 
photograph allowed computation of the scale factor in each case. 

Both light field and dark field photographs were made of the Iso- 
chromatlc fringe patterns in each beam using monochromatic light with a 
wavelength of 5461 A. The resulting negatives were 2.07 times actual 
size, but further enlargement was done so that the total magnification 
of the final prints was about five times actual size in order to facili- 
tate measurement of the fringe spacing. 

A plot of fringe order vs. distance from the neutral axis of the 
beam was made; then a best fitting straight line was drawn through these 
points. The intercepts at the ordinate distances representing the upper 
and lower edges of the beam yielded a total number of fringes correspond- 
ing to twice the fringe order at the extreme fibers. Calculations to 
find f, the material fringe value for stress, were then made. From 
photoelastic theory: [2] 

a:-<r^ . (LL 'i,!' 

where (T, and are principal stresses at a 
point, psi 

N = fringe order at a point 
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f * material fringe value for stress. 



lb/ ir^ringe 

t > model or slice thickness, in. 

From the equation for the extreme fiber stress in a beam in pure bending; 



(T = 



M W 



4.2 



2 I 

where (f = extreme fiber stress 
h = beam height 

I = moment of inertia of the beam cross 
section about the neutral axis 



From 4.2, since I = th "* 

12 



(T = 



6M 



4.3 



This is a tensile (or compressive) principal stress in the outer fibers. 
Since the other principal stress is zero for this case, equation 4.1 
becomes 



or 



(cM 

4 = 



N£ 

t 



4.4 



Since this N is half the value of n, the number of fringes obtained from 



the photograph, the final calibration equation is 

f = 



/2-M , 



Figure 6 shows a typical set of calculations for the material fringe 
value. 

The manufacturer quotes an average value of 1.35 . — — for 

® in-fringe 

Hysol 4290 at a stress freezing temperature of 270°F. One investigator, 

Pih, found values of 1.354 and 1.362 Ib/in/fringe respectively for two 
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Figure 6. TypiCdl Calibration Calculation 






different specimens at a temperature of 280®F. (4] The writer found an 
average value of 1.37 Ib/in— fringe at 280®F. for four specimens from 
the same piece of stock. This value of the material fringe value was 
considered reliable and was used in subsequent analyses of the model 
slices. 

According to Dally, the so-called "critical temperature" for stress 
freezing is not at all critical. [ 2 ] A temperature at or above the 
critical temperature assures rapid response of the material upon applica- 
tion of a load, while temperatures below the critical temperature simply 
result in a time delay in the response. The stress freezing tempera- 
ture for this investigation was chosen to be 280®F. in order to be just 
above the manufacturer's quoted average critical temperature in order to 
assure rapid response of the models irrespective of the uncertainties 
in temperature measurement and variations of the actual critical tempera- 
ture from the average value quoted by the manufacturer. 

The writer made all calibration fringe photographs Immediately upon 
removal from the oven in order to avoid time-edge effects. It was found 
that storage in mineral oil was effective in preventing buildup of time- 
edge effect stresses. 

Section 5. Test Apparatus 

The heating and cooling cycle was accomplished in a Central Scienti- 
fic Company constant temperature oven. Since no time-temperature con- 
trol was available it was decided to use a temperature set point of the 
maximum temperature of 280°F., and to allow the oven to heat and cool 
according to its natural time constant. The empty oven was found to 
heat and cool too rapidly for stress freezing, but the writer found that 
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a container holding about one gallon of mineral oil provided even heating 
of the models and also slowed the heating and cooling rates to acceptable 
values. See Figure 7 for a typical heating and cooling curve for the 
oven-oil tank combination. The oil tank was a stainless steel beaker 
about 8-% Inches in diameter and about 11 Inches high. It easily contain- 
ed the model being tested and a calibration specimen mounted on Its load- 
ing fixture. 

The internal pressure loading of the models was applied by the hydro- 
static pressure output from a dead weight gage tester. This Instrument, 
manufactured by Crosby Steam Valve and Gauge Company, was capable of 
pressures from 0 to 300 pslg in 5 pslg Increnents. The pressure actually 
applied to a model was 10 pslg. This pressure provided a number of 
fringes In the slices large enough to be useful, but the loading was small 
enough to avoid large deformations of the models. 

The connecting piping from the dead weight tester to the model was 
led Into the oven through an existing thermometer hole In the top. 

Figure 8 shows the oven and the general arrangement of the loading piping 
before it had been led into the oven. Since the mercury thermometer was 
removed to make way for the loading piping, a thermocouple junction was 
placed within the central cavity of the model. The thermocouple wires 
were led from the model back through the Inside of the loading piping to 
exit Immediately outside the oven through a small stuffing box. The 
thermocouple system is readily visible in Figure 8. 

A filling and vent connection was placed in the high point of the piping 
in order to fill the system with oil and to bleed off any entrapped air. 

Some slight leakage was present around the thermocouple wires at the 
first application of pressure. After all the air had been expelled, the 
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leaking stopped. Expelling all the air before running a test was essenti- 
al, since the total volume of oil In the loading system was small and the 
distance of travel of the loading screw was limited. With no air present 
In the system the load could be maintained with little movement of the 
loading piston. Otherwise, with air present, the end of the loading piston 
travel was reached before the cooling cycle was complete, and some hur- 
ried maneuvering was required in order to refill the oil reservoir while 
maintaining the required loading pressure on the model. 

Section 6. Test Procedures 

The model was filled with colorless mineral oil, then It was placed 
In the oil tank and connected to the loading piping. The system was fill- 
ed with oil and Initially pressurized to 20 psig to check for tightness. 
This load was removed and the vent plug was removed to allow for expansion 
of the oil during heating; thus, premature loading of a model was avoided. 
The desired maximum temperature for the cycle was set into the oven 
temperature control and heating commenced. After about 9 hours the 
equilibrium temperature, 280° F.,was reached and the test load, 10 pslg, 
was applied. About 45 minutes were allowed for the model response to 
reach equilibrium with the applied load, then the oven was turned off 
and allowed to cool. After the model had cooled, the load was removed 
and the model was ready for slicing and study. 

The manner of slicing a model was determined by considering Its 
geometry and symmetry. Each model consisted of three short cylinders 
having the same axis of symmetry and having junction planes normal to 
the axis. In cylindrical co-ordinates the three normal stresses acting 
upon an element are (T^ , the radial normal stress, , the circumferen- 
tial normal stress, and , the longitudinal normal stress. A thin 
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slice taken in a plane containing the axis of a cylinder contains only 
two principal stresses; at points away from a discontinuity these co- 
incide with <T^ and (Ty. . A slice taken in a plane normal to the axis 
also contains two principal stresses; these coincide with and <T^ . 
Slices taken In a manner to coincide with the plane determined by two 
principal stresses may be analyzed In the same manner as two dimensional 
models, a far simpler process than Is required for a more general three- 
dimensional stress analysis. The slicing plans, therefore, called for 
longitudinal slices and ring slices, each of which would contain a pair 
of principal stresses as described above. 

The slicing plans for the models are shown In Figures 9 and 10. 

Model I had three longitudinal slices removed at 120 degree Intervals as 
shoim In the upper diagram of Figure 9 (slices 1, 3 and 4). These were 
made very nearly the same thickness so that Inspection of the fringe 
patterns afforded an easy check for symmetry of response of the model 
to the applied loading. Figure 11 shows the comparison of these patterns. 
Model I also had three ring slices removed as shown In Figure 9. Model 
II was sliced In similar fashion as shown In Figure 10, but only one 
longitudinal slice and two ring slices were required. 

Slicing was accomplished by using a sharp hand hacksaw with a blade 
having 32 teeth per Inch. Slow careful sawing created only slight dis- 
turbances to the frozen fringe pattern. This distortion was automatical- 
ly removed during subsequent finishing operations. Because of the small 
size of the slices, finishing was easily accomplished by hand sanding, 
using a succession of grit sizes, ending with worn 220 grit aluminum 
oxide paper. Saw marks were removed, coarse scratches eliminated, and 
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the slice thickness was adjusted until It was constant within .001 Inch 
as measured with a micrometer caliper. Polishing was not required since 
an Immersion tank having glass sides was used for containing the slices 
In the field of view of the polariscope. Preliminary tests showed that 
colorless mineral oil having an index of refraction of 1.4775 provided 
a close enough match to the epoxy so that photographs of the fringe pat- 
terns of slices Immersed In this oil would be satisfactory. Figure 12 
shows the Immersion tank being used with the polariscope. 

The polariscope, Figure 13, Is a Chapman condensing type having a 
5 Inch field and an Integral camera with an 8x10 Inch back. Two light 
sources are Interchangeable: one white light source and a mercury vapor 

source filtered to provide monochromatic light of 5461 A wavelength. 

The camera arrangement makes negative Images that are 2.07 times the 
actual size of the object. Contact prints of the fringe patterns were 
therefore of convenient size without additional enlargement. 

Kodak Contrast Process panchromatic film was selected because of 
Its high contrast and Its good sensitivity to green light. Exposure 
times were about 1 to 2 seconds with wide open aperture. Fringe photo- 
graphy was done with open aperture because it provided an evenly illumin- 
ated field. Small apertures tended to produce a partial image of the 
light source superimposed upon the fringe pattern; this effect produced 
poor negatives. In general, both light field and dark field photographs 
were made of the model slices and calibration specimens. 

Fringe orders were most conveniently found by the method described 
by Dally. [2] A wedge was cut across any convenient portion of the 
slice. A zero order fringe then appeared on the thin end of the wedge. 
This method was particularly applicable to the ring slices, where all 
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fringes were of high order. For the longitudinal slices, a free exterior 
comer provided a zero order starting point. A third method was to 
employ a stress frozen calibration beam as a compensator. This last 
method was found to be a useful supplement to the other two. 

Section 7. Results of Tests 

The objective of the investigation was to study the stress distri- 
butions for the two geometries in order to determine what differences 
existed as a result of the variations in geometry. 

The simple body-of-revolution geometry allowed correlation with the 
Lame^ solution for elastic stresses in thick -walled cylinders with internal 
pressure loading. The derivation of this solution is available from many 
texts; for example, see Dally. [2] The results of this solution are equa- 
tions for (T^ and , the radial and circumferential principal stresses, 
respectively. Figure 14 shows the diagram for the solution and indicates 
the meaning of the symbols employed. 

The circumferential principal stress is: 




The radial principal stress is: 




where = internal loading pressure, psig 
a = inside radius 
b « outside radius 
r =* radius (variable) 

Equations 7.1 and 7.2 may be combined to yield a stress difference equation: 
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Figixre 1^, Diagram for The Lame^ Solution For Principal 
Stresses In A Thick Walled Cylinder. 
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Note that the photoelastic equation for a ring slice takes the form 









7.4 



t 



where N = fringe order 



f = material fringe value for stress, Ib/in/fringe 



t = slice thickness, in. 



Now the experimental values for the stress difference could be compared 
with the theoretical values in order to check the reliability of the 
test results. The results of this comparison for the ring slices from 
Model I are shown in Figures 15, 16, and 17. The experimental points 
agree rather well with the predictions of the theory. 

It should be noted that equations 7.3 and 7.4 may be combined and 
rearranged; 



This relation may be used as a quick check on the consistency of a fringe 
photograph of ring slices, for the right hand side is a constant for a 
given slice. 

Unfortunately, no such simple result was available to use for check- 
ing the longitudinal stress, (f^ , or the longitudinal stress difference, 

<T^ • As mentioned in Section 2, simple cylindrical shell theory 

was not applicable because a third normal stress, (3^ , was present in 
addition to the two membrane stresses, and . This radial stress 
was not negligibly small; in fact it frequently had a magnitude comparable 
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Figiire 16 . Model I, Slice Comparison of Experimental 
and Theoretical Values of C-c. 
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Figure 17. Model I, Slice 6. Comparison of Experimental 
and Theoretical Values of 6L-Ti-» 



to that of (Tg^ . Another factor causing difficulty was the presence of 
an edge moment at the junction of the smaller, 1.900 in. diameter cylinder 
and the larger, 2.200 in. diameter cylinder. This moment was the result 
of the change in the mean diameter of the cylinders at the plane of 
junction. Another effect resulted from the differing magnitudes of the 
radial and circumferential strains in each of the two cylinders at points 
just away from the junction plane. Compatibility requires that these 
strains be equal in the junction plane; this can be the case only with the 
existence of additional shearing forces and edge moments within the cylin- 
ders at the discontinuity. 

A solution could be devised to this complicated problem by using 
some finite element numerical analysis scheme programmed for a digital 
computer, but this method of analysis is beyond the scope of this paper. 

The good agreement in Model I of the experimental values of the 
stress differences in the ring slices with those obtained from the Lame 
solution suggested that the theoretical value of <7^ in each case could 
be used as a method of separating the stresses in the longitudinal slices 
at locations where and were the principal stresses. Determina- 
tion of by itself would make it possible to obtain a static check. 

The procedure followed in this process was first to obtain a fringe 
order vs. radius curve. Multiplying the order by the slice fringe value 
furnished the value of the stress difference at each point. The Lame' 
Bolution curve for was plotted on the same axes. The difference in 

ordinates between the (Q-(T^ curve and the curve was plotted as the 
curve for C7^ . This last result was then used to plot a curve of 

vs. r. The static check at a section required a summation of the normal 
stress over the area of the section, an integration of the form 
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2-n 




(T^r dr 



7.6 



\^ere P ■ “ the normal force at a section. The value of the 

Integral In equation 7.6 was found graphically by applying Simpson's 
rule to find the area under the (T^f vs. r curve. This area was multi- 
plied by 2^ and the graph scale factor for the area units to yield the 
value of P. The appropriate graphs and the calculations for Model I 
may be found In Appendix II. For statics to be satisfied this value of 
P had to be Identical to the opposing force, the result of the Internal 
pressure acting on the Inside cross sectional area of the cylinder. The 
result of this check for Model I Is shown In the following table. 



SECTION 


STATIC CHECK LOAD. Lb. 


APPLIED LOAD. Lb. 


% DIFFERENCE 


Pine 


18.3 


17.3 


+ 6 


Coupling 


16.3 


17.3 


- 6 



The differences are considered acceptable for a static check. 



The curves of 0^vs. radius were useful for determining the behavior 
of the model, but with, certain reservations. The stress distribution 
generally followed what was expected considering that shell bending was 
superimposed upon the normal Lame^ stresses. For example. In the 
curve for the zone of slice 5 (the pipe) In Figure 25, ^ Is seen to be 

largest In the outer fibers and smallest Inside. This follows, for the 
shell bending moment at this section must have produced tensile stresses 
In the outer fibers and compressive stresses In the inner fibers, all 
superimposed on a constant (7^ of about 16.3 psl. Similarly, In the 
coupling section at slice 6 the ^ curve shown In Figure 28 shows some 
bending effects In the opposite sense, for the largest values are near 
the Inner surface of the coupling. In these curves the photoelastlcally 
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determined values of <f^ were low at the boundaries. These values are 
probably erroneous for two reasons: first, it is possible that the "time- 

edge effect" superimposed local compressive stresses upon the already ex- 
isting tensile system; or second, that the low indicated fringe order near 
a boundary may have resulted from the variation in slice thickness at the 
boundary, a natural consequence of the fact that the surfaces at the 
boundaries are portions of cylindrical surfaces rather than planes. If 
these edge effects are ignored, the other indicated values of (T^ appear 
to be believable. 

The longitudinal slices, numbers 1, 3, and 4, showed in their fringe 
patterns that there was some sort of stress concentration taking place at 
the discontinuity of the junction, but the slice geometry and character- 
istics of the cement produced rather broad shadows which effectively 
obliterated the fine pattern which may have existed at the junction of 
the cylinders. For this reason no stress concentration factor was obtain- 
ed for this discontinuity. 

The slight difference in the ends of the coupling proved to be in- 
significant. It did not alter the stress distribution except to make a 
somewhat smoother transition between the small cylinder and the large 
cylinder at the end having the 45 degree bevel. 

The 'analysis of Model II was conducted in similar fashion, but the 
results were erratic. Figures 35 and 36 in Appendix 111 show the compari- 
son of the photoelastic stress difference and the theoretical Lame^ 
curve for the coupling and pipe ring slices. For the coupling slice 
the experimental values are all low by what appears to be a nearly con- 
stant amount. The single experimental point for the pipe slice fell 
remarkably close to the theoretical curve, making the two sets of results 
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mutually inconsistent. Whatever affected slice 2 apparently did not 
affect slice 3. Another check on the experimental results was to attempt 
to satisfy statics as was done for Model I. The ring slice results did 
not completely justify the assumption that the Lame'' values for (T^. could 
be used to separate the principal stresses, but a static check was made 
using this assumption as an approximation. The stress distribution for 
the section at the center of the coupling is shown in Figure 38, and the 

r vs. r curve with Its accompanying force calculations are shown in 
Figure 39 and following In Appendix III. The calculated force here was 
only 11.6 lb., 33% lower than the applied load of 17.3 lb. This tended 
to confirm that the results of Model II were too much In error to be 
quantitatively useful. 

Although good numerical results were not obtainable from the test 
of Model II, some Inferences may be gathered from Model I about the In- 
fluence of the circumferential groove on the stress distribution. The 
fringe pattern of slice 1 from Model I (Figure 19) shows the distortion 
of the pattern at the groove that corresponds to the Increase In level of 
the stress intensity at the narrowing section. It appears that extend- 
ing the fringe pattern on either side of the groove would yield a pattern 
corresponding to the coupling with no groove. Alternatively, the results 
of the analysis for the zone of slice 6 could be used to represent the 
plain coupling. Figure 40 shows the variation of (7^ with radius for the 
section taken at the groove. It shows an approximate mean stress of 15 
psl, while the mean stress of the undlmlnlshed section is about 9 psl. 

The ratio of these stresses Is about 1.7. The corresponding Inverse ratio 
of the cross section areas Is about 1.4. This shows very roughly that the 
stress level In the section at the groove Is affected largely by the decrease 
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in section area. 



Section 8. Suggestions for Further Investigation 

A cemented model such as was used to represent brazed joints is not 
very good for studying the effects of discontinuities. A one-piece model 
with cast or cemented ends would be better because the effects would be 
much clearer. Cement lines, shadows and attendant air bubbles would all 
be eliminated. The geometry at comers and fillets would be easier to 
control because of the absence of cement in these important areas. 

Determination of boundary stresses depends upon freedom from the 
time-edge effect. In the tests described, even very rapid handling of 
the models and slices was not adequate to avoid some interference from 
this effect. 



43 



BIBLIOGRAPHY 



1. Burke, RoberC E. Photoelastic Study of Pressure Vessel Stresses . 

Unpub. M. S. Thesis, U. S. Naval Postgraduate School, 1955. 

2. Dally, J. W. and Riley, W. F. Experimental Stress Analysis . 

McGraw-Hill, 1965. 

3. Hetenyl, M. Ed. Handbook of Experimental Stress Analysis . 

Wiley, 1950. 

4. Plh, H. "Three Dimensional Photoelastic Investigation of Circular 

Cylinders with Spherical Cavities In Axial Loading." 
Experimental Mechanics , v. 5, March, 1965; 88-95. 

5. Timoshenko, S. and Wolnowsky-Krleger, S. Theory of Plates and 

Shells, 2nd Ed . McGraw-Hill, 1959. 



44 



APPENDIX I 



MATERIAL AND MODEL DATA 
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Technic ctl Data for Hysol 4290 
Typical Properties 



Critical temperature 

Tensilie strength @ 77° P 

(§ 270° F 

Modulus of elasticity @ 

@ 270 F 

Poisson's ratio - 

Fringe value @ 77° ^ 

@ 270° P 

Annealing temperature 

Model I Actual Dimensions 
Inside- diameter 
Pipe outside diameter 
Coupling outside diameter 



270® F 

12000 psi 
210 psi 

480000 psi 
2190 psi 

,500 

57 Ib/ln-fringe 
1.35 Ib/in-fringe 

280'’ - 300° P 



1.500 in. 
1.898 in. 
2.177 in 



Model II Actual Dimensions 
Inside diameter 
Pipe outside diameter 
Coupling outside diameter 



1 . 500 in* 

1.891 in. 

2.177 in. 
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Calibration Beam Moment Calculations 



Applied Loc*d: 


lb. ' R- 

|er 4<^o' >j 


: ,175 lb. 








' 1 



•ns li. 



Load Diagram 



.ns lb. 



+. n^lh 
a 



-./7r lb. 



Shear Dia,^ram 



0 






Moment Dia«;ram 



Gravity Load: 



Density of epoxy 
Density of oil • 
Apparent weight of beam 
Load 



.053 Ib/in^ 
.030 Ih/ln^ 
.021 lb 
.002^ Ib/in 



Hr 



,OIOZ lit 






. »» 
ici> 









..V y y _V_v .V 



y V \:/ 

tv 






o 



fb/tn. 



Load Diagram 

f,Oo^U> 




Shear Diagram 



,00^7 

I-bment Dla-^rain 




0 



^, 00 '^y 



(Maximum moment in the center span is less thcin 
\$ of the applied moment.) 
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APPENDIX II 

MODEL I FRINGE PHOTOGRAPHS AND CALCULATIONS 



■ t ' 
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I 




I 

I 4^ 

■ ■ I 
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Fig jre 23- 



' odcl I, oiioc Fi'l'ige P.-ittor' s for Dar!; 
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MODEL I calculations FROM PHOTOGRAPH DATA 



Slice 2 






fc lf87^7.2 psi/fringe 




N -r,in. 


^ ^ ^ ^ psi 


(Lame') ,psi 


3.5 *982 


25.2 


26.6 


^ .938 


28.8 


29.2 


^.5 .873 


32.4 


33.6 


5 .806 


36.0 


39.5 


5.5 .777 


39.6 


h2.7 


Slice 5 



T S' ljt3Z-9.62 psi/fringe 
T TIlf2 



N r, in. 


psi 


(Lame) , psi 


4 .926 


38.4 ^ 


34.8 


4,5 .866 


43.3 


4o.o 


5 .784 


48.1 


48.6 


Slice 6 



f «li3Z-7.8 psi/fringe. 
t '7175 



N 


r, in. 


< 5 ; , psi 


<Tg--(r^ (LameO, psi 


2 


1.04 


15.6 


19.9 


2.5 


.992 


19.5 


21.9 


3 


.938 


23.4 


24.4 


3.5 


.902 


27.3 


26.4 


4 


.844 


31.2 


30.5 


4.5 


.798 


35.1 


33.8 


5 


.775 


39.0 


35.8 



\ 
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MODEL I LONGITUDINAL STRESS CALCULATIONS 
FROM PHOTOGRAPH DATA 



Slice 4 



,7.2 psi/frince 



N 


r, in. 


, psi 


psi 


psi 


‘^z^,lt>/in 


2 


.750 


l4.4 


-10.0 


4.4 


3.3 


2.5 


.786 


18.0 


-7.6 


10.4 


8.2 


3 


.850 


21.6 


••4.0 


17.6 


14,9 


3.5 


.910 


25.2 


-1.4 


23.8 


21.6 


3 


. 9^9 


21.6 . 


0 


21.6 


20.5 


^ne 


6 










N 


r, in. 


psi 


psi 


^4 psi 


Ib/in 


2.5 


.750 


18.0 


10.0 


8.0 


6.0 


3 


.788 


21.6 


8.2 


13.4 


10.5 


2.5 


.832 


18.0 


6.4 


11.6 


9.6 


2 


.890 


l4.4 


4.5 


9.9 


8.8 


1.5 




10.8 


2.8 


8.0 


7.6 


1 


1.028 


7.2 


1.0 


6.2 


6.4 


tl— 


1.064 


3.6 


0.4 


3.2 


3.4 


Zone 


_2 _ 










N ■' 


r, in. 


psi 









2.5 ,/96 18 

3 .77 21.6 
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Figure 2^. Ilodel I, Slice V. Fringe Order vs. Radius at Zone 
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Figure 25, ^-lodel I, Slice h. Stress Distribution at Zone 




I 



$ 

o 



\ 



58 



Figlire 26, Model I, Slice h,’' Curve for Graphical Force Sumation 



MODEL I STATIC CHECKS 



Zone 5 (See Figure 26.) 

Simpson’s rule was applied to the O^'C vs. r 



The abscissa was divided into 8 segments of 1 inch 
on the graph, having a scale value of .025 inch. 
This yields an area too large by the amount of the 
excess included in the last segment; call this area 



a' St 1 j”(,664- 4.1) + 4(1. 4+2. 4+3. 7f 4.4)42(1.9+2.9+4.2) 



A' =s 23.53 sq. in. 

A 23.37 sq. in. 

Graph area value - .125 lb. 

/, P=(27I)(23.37)(.125)s:l8.3 lb. 

Considering internal pressure loading on inside cross 
section: P — ^a^ 

P=(10)Cr\)(.75)2 =17.3 lb. 

Per cent difference: 



curve in order to find the value 




a 1 Then the required area is A= A'- .16 square inches. 



18.3-17.3 X 100=-+6^ 

17.3 











6Q 



Figure 27 • 1-odel I, Slice Fringe Order vs. Radius at Zone- 
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Figure 28, Model I, Slice Stress Distribution at Zone 
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Figure 29. Model I, Slice 4. Curve for Graphical Force Sumracition 
at Zone 6, 



MODEL I STATIC CHECKS 



Zone 6 (See Figure 29 •) 

Simpson's rule WdS applied with the abscissa divided 
into segments 1/2 inch long. The area unit used was 
a square 1/2 x 1/2 inch having a scale value of .0625 lb. 
A** A — .08 sq. units 

a'= ir(2. 44. 7+3. 043. 6+3. 3-*-2. 9+2. 641.0) , : 

3L 1 

+ 2(4.2+3.8+-3.5+3.042.8+1.8)J 

a^-4i. 67 sq. units 
A -41.59 sq. units 
Graph area value =-.062 5 lb. 

P= (4i. 59)(27T)(.0625) =^16.3 lb. 

Considering internal pressure loading on the inside 
cross section: P=r- p./ra^ 

P» (10)(7r)(. 75)^= 17.3 lb 

Per cent difference; 

16.3 -17.3 X 100- -6^ 

17.3 
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appendix III 

MODEL II FRINGE PHOTOGRAPHS aND CALCULATIONS 
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TG 30* 2'odol II. Slice 1. Frin.'^e Patterns for Dari: and Light Field 




3^* ocIgI II5 0Z.J.CG 



o. Tri t G c tto:- fo 2 



T*' * r"! 
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MODEL II Calculations from photograph data 



Slice 2 



f — 1 « ~^7 — 6.06 psi/fringe 

t .226 



N 


r, in. 




, psi (T^-Tr (IfcinieO 


, psi 


2.5 


1.03 


15.1 


20.3 




3 


.95 


18.2 


23.9 




3.5 


.90 


21.2 


26.6 




If 


.84 


24.2 


30.5 




5 • 


.80 


27.2 


33.6 




CO 

o 

CD 


3 








f 




6.79 p si/fringe 






t 


.202 






N 


r, in. 




psi - (T^(Lajne'" 


), psi 


6 


.85 


40.7 


41.4 


> 


Slice 


1, zone of 


slice 2 






N 


r, in. 


, psi 


( 

^9 psi d^, psi 


<^r, Ib/in 


1 


1.04 


6.85 


.9 5.95 

3.6 6.7' 


6.2 • 


1.5 


.92 


10.3 


6.2 


2 


.78 


13.7 


8.5 5.2 


4.1 



I 
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Figure 35 * I-odel II, Slice 2. Comparison of Sxperiraenty.l and 
Theoretical Values ot iJq ~ , , 
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Model II, Slice 1. IJi*inge Order vs. Radius at Zone 
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Figure 38. Model II, Slice 1 . Stress Distribution at Zone 




‘o 

I 



§ 



CM 



75 



Figure 39. Model II, Slice 1. Ctu^ve for Graphica.! Suinraation at Zone 



MODEL II static CHECK 
Zone 2 (Figure 39) 

Simpson's rule was applied with the abscissa divided 
into l4 segments 1/2 inch long. The area unit used was 
a square 1/2 x 1/2 inch having a scale value of .02? lb. 
As A — 2.48 sq. units 

i r (3. 4-<6.2H4(4. 0+5.0^5.846. 246. 2+6. 2+6. 2) 

+ 2(4. 5+5*446. 1+6. 2+6. 2+6. 2)/ 

A ^=76.4 sq. units ,, ' • 

A ■= 73.9 sq. units 
P = (73.9)(2jT)(.025)^ 11.6 lb. 

Considering internal pressiire loading on the inside 
cross section area: P = p^na^ 

P=(10)(7T)( .75)2 . ,7,3 
Per cent difference: ^ 

11.6-17.3 X 100 =-33^ ' 

17.3 
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appendix IV 





Figiire 40. Model I, Slice 1. Stress Distribution at the Coupling 
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